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By analytical mapping of the eigenvalue problem in rough billiards on to a band random matrix 
model a new regime of Wigner ergodicity is found. There the eigenstates are extended over the 
whole energy surface but have a strongly peaked structure. The results of numerical simulations 
and implications for level statistics are also discussed. 

PACS numbers: 05.45.+b, 72.15.Rn, 03.65.Sq 



OS 
PL, 



> 

in 

CO 

CN 

O 

r^ 

o> 

-i— > 



i 

C 

o 
o 



* 



In 1974, Shnirelman |l| proved a theorem according to 
which quantum eigenstates in chaotic billiards become 
crgodic for sufficiently high level numbers. Later it was 
demonstrated [g,f3| that in this regime the level spacing 
statistics p(s) is well described by random matrix the- 
ory H- However, one can ask the question how this 
quantum ergodicity emerges with increasing level num- 
ber iV? This question becomes especially important in 
the light of recent results ymw for diffusive billiards where 
the time of classical ergodicity td due to diffusion on the 
energy surface is much larger than the collision time with 
the boundary t&. In such a situation quantum localiza- 
tion on the energy surface may break classical ergodicity 
eliminating the level repulsion in p(s). The investiga- 
tion of rough billiards |6| showed that this change of p(s) 
happens when the localization length I in the angular 
momentum /-space becomes smaller than the size of the 
energy surface characterized by the maximal I = l max at 
given energy (I < l max ). For t > l max the eigenf unctions 
are extended over the whole surface but as we will see 
they are not necessarily ergodic (Fig. [j]). 

In this situation which we will furthermore call Wigner 
ergodicity the eigenstates are composed of rare strong 
peaks distributed on the whole energy surface. Such 
a case is very different from the Shnirelman ergodicity 
where the eigenstates are uniformly distributed. The 
usual scenario of ergodicity breaking was based on the 
image of quantum localization along the energy surface 
J7J. Here we show that the transition between localized 
and Shnirelman ergodic states can pass trough an inter- 
mediate phase of Wigner ergodicity. Our description and 
understanding of this phase is based on the mapping of 
the billiard problem with weakly rough (random) bound- 
ary on to a superimposed band random matrix (SBRM). 
This model is characterized by strongly fluctuating di- 
agonal elements corresponding to a preferential basis of 
the unperturbed problem. Recently such type of matri- 
ces was studied in the context of the problem of particle 
interaction in disordered systems M flj. There it was 
found that the eigenstates can be extended over the whole 
matrix size while having a very peaked structure. The 
origin of this behavior is due to the Breit- Wigner form 
Jl2[ of the local density of states according to which only 



unperturbed states in a small energy interval Te con- 
tribute to the final eigenstate. 




FIG. 1. Transition from localization to Shnirelman ergod- 
icity on energy surface for level number N ~ 2250, l ma x ~ 95 
and M — 20; shown are the absolute amplitudes |C„j | of one 
eigenstate: (a) localization for D(l r = 0) = 20; (b) Wigner 
ergodicity for D = 80; (c) Shnirelman ergodicity for D = 1000 
(see text). 



Recent optical experiments with micrometer size 
droplets initiated new theoretical investigations of weakly 
deformed circular billiards |13|. In this case the ray dy- 
namics becomes chaotic leading to a strong directional- 
ity of light emission |L4| . Here we will consider another 



type of a weakly deformed circle Q, namely we chose 
a random elastic boundary deformation which can be 
represented by R(6) = R + AR(0) with AR(6)/R = 



,1/ 



with Bessel functions Ji and i? being a normalization 
constant. The S- matrix gives a local unitary description 
for I close to a resonant value Z r . The eigenvalue equation 
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where 7 m are random complex co- reads Y^i^iii-^a) 
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efhcients and M is large but finite. This type of de- 
formation seems to be very generic and may appear 
in numerous different physical situations |6|. In small 
droplets such boundary perturbations may be created 
by temperature induced surface waves. We will restrict 
ourself to the case of weak surface roughness given by 
k(6) — (dR/d8)/R <C 1 and all -f m being of the same 
order of magnitude. Then we have for the angle average 
k 2 = (n 2 {9)) 9 ~ M 2 (AR/R ) 2 . 

In such a billiard the dynamics is diffusive in or- 
bital momentum due to collisions with the rough bound- 
ary provided k is above the chaos border k c ~ M~ 5 ' 2 
P|. The diffusion constant is determined by the av- 
erage change of orbital momentum per collision being 
D =< (A/) 2 >= 4 {l 2 nax -l 2 r ) k 2 . This D is the local diffu- 
sion rate for I close to l r . The quantum interference leads 
to localization of this diffusion with the length £ = D for 
M < £ < l max while for I > l max the eigenstates are 
extended over the energy surface |]6|. The transition be- 
tween these two regimes is illustrated in Fig. [y. Here 
we present the absolute values of eigenfunction ampli- 
tudes C„" in the cigenbasis \nl > of circular billiard as a 
function of unperturbed radial and orbital quantum num- 
bers n, I with a marking the eigenenergy E a . For small 
roughness k (or D) the states are exponentially local- 
ized (Fig. |lja) while for large k they are homogeneously 
distributed (Fig. [j]c) on the energy surface. The case 
of Fig. |ljb corresponds to an unusual regime of Wigner 
ergodicity where the eigenstate is extended over the sur- 
face but is composed of rare strong peaks. The positions 
of these peaks on the energy surface of the circular bil- 
liard E — TL(n, 1) are shown in Fig. ga. The equation of 
the surface, projected on the action plane (n,l), can be 
found from the Bohr-Sommcrfcld quantization fJ-i(E) = 
V l la X ~ I 2 ~ I arctan(r x ^Jl 2 max - I 2 ) + tt/4 = n(n + 1) 
where l 2 nax = AN = 2mR 2 E/h 2 = k 2 R 2 with k being 
the wave number. A part of the surface is shown in more 
detail in Fig. ||b. Here it is clearly seen that the peaks 
are large for those integer n, I which are close to the 
line Ti(n,l) — E a . Our understanding of the fact that 
not all integer values of the (n, Z)-lattice near to this line 
are populated is based on the concept of Breit- Wigner 
structure of eigenstates described below. 

According to Refs. IqJlql the internal scattering at the 
rough boundary can be described by the S'-matrix 



S d (E) = e Hll(E) < l\e lV{s) \l> e Hl i {E) , 



(1) 



with Hi (E) being the scattering phases of the circle and 

V(0) 

map 



2 V l2 nax ~ l r AR{9)/Ro. This quantum rough 
is defined with respect to amplitudes a; in the 
wave function expansion ip{r,6) = B^2 t ai Jm(kr) e lW 



so that the eigenvalues E a 
are determined by det[l — S(E a )] — 0. For V = 0, we re- 
cover the Bohr-Sommerfeld quantization for eigenvalues 
E n i of the ideal circle. 
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FIG. 2. (a) Main peaks of eigenstate in Fig. lb (squares 
> 0.1) shown on the energy surface TC(n,l) = E a 
(b) rescaled part of (a): diamonds show the integer 
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(n, Z)-lattice, the errorbar size is 2|C, 

The semiclassical regime of ray dynamics corresponds 
to the limit V 3> 1 where the ^-integral can be evalu- 
ated in a saddle point approximation giving the classi- 
cal limit of quantum rough map 0. Here we are inter- 
ested in a different regime where V < 1 corresponding 
to D < At 2 . There, by the mapping on an effective solid 
state Hamiltonian H e ff introduced by Fishman, Grcm- 
pel, and Prange Q, the equation for eigenstates takes 
the form 



ta,n\jj,i(E a )]a 



(a) 
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(2) 



In this way the eigenvalue equation is reduced to a solid 
state problem with 2M coupled sites. The if e jy-matrix 
is of SBRM type with strongly fluctuating diagonal ele- 
ments produced by scattering phases /z/. The investiga- 
tions of such matrices [p|-pT[ showed that the local density 
of states has the Breit- Wigner width given by the Fermi 



golden rule r A 



2w Pli < (V{6)/2) 2 >« 3L»/(2Af 2 ) 



where p M = 1/-7T is the density of diagonal elements and 



we used the relation between phase average of V 2 (9) and 
D. This expression is valid |17j when T M exceeds the 
mean level spacing (~ 1/M) in the band width M. In 
the opposite limit T^M < 1 the eigenstates are given 
by standard perturbation theory. Together with the con- 
dition V < 1 we find that the Breit-Wigner regime ex- 
ists for M < D < M 2 near zero energy of Heff- In 
this regime the localization length is £ = D [||,|| |Ii||l5| |. 
However, the Breit-Wigner structure remains in both lo- 
calized (£ < Irnax) and delocalized (£ > l max ) cases if 
M < D < M 2 . Therefore for l max < D < M 2 the states 
are extended but only / with |tan(/^)| < T^ < 1 are 
mixed leading to a peaked structure of eigenstates [[18| . 
The fraction of peaks in max(l, lmax) is T M . 

The above properties of scattering amplitudes a\ al- 
low to understand the behavior of eigenfunction coeffi- 
cients C^ —< ipa\nl >. For this one has to compute the 
expansion Ji(k a r)e lW in terms of \nl >. Since AR <C Ra 
the angular and radial integrals factorize and can be eval- 
uated using the radial eigenvalue equation and the semi- 
classical expression for Ji(kr). As a result we obtain 
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with A/i = fii(E a ) - [ii(E n i) sw ( Eg - E n i )/E b and 
E b = dE/dvLi{E) = h 2 l 2 ma J(mR 2 ^ll ax -l 2 ) = 2h/r b 
being the energy scale related to the ballistic collision 
time Tb, B is a normalization constant. The amplitudes 
C^" determine the local density of states by 



pw(E - E n i) = (J2 S ( E ~ E ^ \ C : 



(«)|2 

I I 



(4) 



The averaging is performed with respect to different 
roughness realizations and/or over a sufficiently large en- 
ergy interval. Due to the Breit-Wigner distribution for 
tan(/i() in (g) we obtain 



Pbw(E — E n i) — — 



T E /2 



it {E-E nl ) 2 +r 2 E /4 



with 



r_E = Eb r„ = Eb 



N 
Nw~ 
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■max 



, N w 



M 2 
24k 2 



(5) 



(6) 



The equations (§), (|) are valid for AE < E B (A/i < I) 
and Tft < 1 or N < Nw- Here Nw is the border of 
Breit-Wigner regime in level number N. We remind that 
the eigenstates are localized for N < N e = 1/(64k 4 ) cor- 
responding to I < l m ax fj) ■ As a result the Breit-Wigner 
structure can exist both in the localized and delocalized 
cases. An example of Breit-Wigner distribution is shown 
in Fig. H. Our numerical data confirm the theoretical 
expression (0) for variation of T^ by more than one order 
of magnitude (inset). 
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FIG. 3. Breit-Wigner distribution for eigenstates of rough 
billiard (diamonds) with the parameters of Fig. lb (5 eigen- 
states for each of 10 roughness realizations are used). Full 
curve gives the distribution (pi) with the theoretical Fe value 
(rj), Vs/Eh — T M = 0.3. Insert shows the variation of 
Ts/Eb (diamonds) as a function of r M for parameter range 
10 < M < 40 and M < D < M 2 (E b and T M are taken at 
I — 0). The theory (@) is shown by straight line. 

For N > Nw the kick amplitude V in ([l]) is larger 
than 1 and the mapping on to equation (0) is not valid. 
In this case the scattering phases (eigenphases of S) are 
homogeneously distributed in the interval (0, 2ir). If in 
addition N > N e then as in the case of kicked rotator 
(sec B. V. Chirikov in 0]) the amplitudes ai are homoge- 
neous in ^-space with \ai\ 2 sw l/(2l max ). Using Eq. (||), 
we obtain that the local density of states is given by 



pw(E - E n i) 



E b sin 2 [{E-E nl )/E b ] 
it [E-E nl f 



(7) 



This density is normalized to one and as a result 
the probability |C„" | 2 is ergodically distributed along 
the energy surface shown in Fig. 0a. This is the 
regime of Shnirelman ergodicity which emerges for N > 
m&x(Nw, N e ). For fixed roughness k > kew — \/6/4M 
we have Nw > N e and the transition to Shnirelman er- 
godicity with the increasing level number N crosses the 
region of Wigner ergodicity for which an eigenfunction is 
ergodic only inside the Breit-Wigner width Ye < Eb- In 
the opposite case k c < k < kew the Shnirelman ergodic- 
ity emerges directly from the localized phase (the Breit- 
Wigner regime exists only in the localized phase). The 
averaging of equation (In) over different Z-values gives the 
distribution (H) with Te ~ Eb. This explains why also 
for the case r M < 1 in Fig. g the distribution (J5|) remains 
valid even for AE > Eb (note that nEb/Vs « 10). 

The above analysis shows that in the regime of Wigner 
ergodicity there are four relevant energy scales: level 
spacing A, Thouless energy for diffusion in Z-space E c = 
hD/(l max T b ), the Breit-Wigner width T E = 3hD/(M 2 r b ) 
and bouncing energy E b = 2h/T b which are ordered as 
A < E c < Te < Eb- These scales should appear 
in the level statistics namely for the number variance 



T, 2 (E) [gjj]. For E < E c as usual we expect gaus- 
sian orthogonal ensemble (GOE) statistics to be valid 
while in the interval E c < E < Te/2 the behavior 
should be modified, due to the diffusive dynamics p9[ , 
being ^(E) ~ (E/Ec) 1 / 2 . The first investigations of the 
regime with T E /2 < E < E b /2 for SBRM were done only 
recently J20| . They showed that level rigidity is strongly 
suppressed with a nearly linear energy behavior in S2 (E) 
due to disappearance of correlations between levels with 
energy differences larger than Te- However, such lo- 
cal characteristics as p(s) are still described by GOE if 
E c 3> A. Our numerical data qualitatively confirm these 
expectations (see Fig. |4|), however, quantitative numeri- 
cal and analytical verifications are still required. In Fig. 
[| the above energy scales are not separated by strong in- 
equalities but parametrically it is possible to have them. 
In this unusual regime it would be interesting to study 
other physical properties. We mention for example the 
frequency dependence of dielectrical response [|2lf which 
should be sensitive to the above energy scales. 



and ip a has a strongly peaked structure on (n, Z)-lattice. 
However, at the same time all states in the Breit-Wigner 
width r^ are populated so that the eigenfunction is er- 
godic inside this energy band. It would be interesting 
to understand if the arithmetical properties of the lattice 
will play an important role. 




FIG. 4. Dependence of number variance E2 (E) on energy 
for rough billiard compared to Poisson and GOE for M = 50, 
D = 800, Imax ~ 95; 6 roughness realizations in the interval 
2150 < N < 2350 are used for the average. The energy scales 
are also shown in units of level spacing A. 

In conclusion, we studied the parameter dependence 
of the quantum energy surface width in rough billiards. 
In the limiting case of Shnirelman ergodicity with high 
level numbers, this width is determined by the typical 
frequency of collisions with the boundary (r^ ~ Ef,). 
This means that all integer points on the (n, Z)-lattice of 
quantum numbers with a distance Al = An w 1 from 
the energy line E a = Ti(n, I) are occupied by one eigen- 
function tp a (N > Nw and N > N e ). We have found a 
new regime of Wigner ergodicity where Te <C Eb so that 
only points with Al — An < T p < 1 contribute to ift a . 
As a result a lot of holes appear in the energy surface 
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